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Abstract. We introduce the extension of the ab initio no-core shell model with contin¬ 
uum to describe three-body cluster systems. We present results for the ground state of 
®He and show improvements with respect to the description obtained within the no-core 
shell model and the no-core shell model/resonating group methods. 


1 Introduction 

The ab initio no-core shell model/resonating group method (NCSM/RGM) was presented in Qi 
as a technique that is able to describe both structure and reactions in light nuclear systems. Within 
this approach, the wave function is expanded in a continuum cluster basis using the resonating group 
method with realistic interactions and a consistent ab initio description of the nucleon clusters. 

The method was first introduced in detail for two-body cluster bases and has been shown to work 
efficiently in different systems iMl . Later, the expansion of the method for three-cluster systems was 
introduced in Si. 

The capability of ab initio methods to properly describe three-body cluster states is essential for 
the study of nuclear systems that present such configuration. This type of systems appear, e.g, in 
structure problems of two-nucleon halo nuclei such as ®He and "Li, resonant systems such as or 
transfer reactions with three fragments in their final states such as ^H(^H,2n)"^He or ^He(^He,2p)"^He. 

Despite the success of the NCSM/RGM in describing the long range behavior of the wave func¬ 
tions, it has been shown that it has limitations when it comes to accurately account for short range 
correlations, which is necessary to achieve a complete description of the system. This is due to fact 
that to account for such correlations, several excited states of the nuclear clusters must be included 
in the basis, resulting in an increase of the problem size that goes beyond current computational ca¬ 
pabilities. This limitation has been overcome by introducing the ab initio no-core shell model with 
continuum (NCSMC). With this method, the wave function is written as a superposition of both con¬ 
tinuum NCSM/RGM cluster states and discrete eigenstates of the compound system obtained with 
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Table 1. Energy (in MeV) for the ^He ground state using the NCSM/RGM, NCSM and NCSMC approaches at 
A^maA-=12. For the NCSM we also show the extrapolated value to N^ax —» o°- 


^max 

NCSM/RGM 

NCSM 

NCSMC 

8 

-28.62 

-28.95 

-29.69 

10 

-28.72 

-29.45 

-29.86 

12 

-28.70 

-29.66 

-29.86 

Extrapolation — 

-29.84(4) 

— 


the no-core shell model (NCSM). The latter eigenstates compensate for the missing cluster excita¬ 
tions improving the description of short range correlations. 

The NCSMC was first introduced in ITl [st] for binary systems. Its expansion to three-cluster 
systems was recently achieved and we show here the first results for the ®He ground state (g.s). 


2 Formalism 


In the NCSMC, the ansatz for the three-cluster many-body wave function is given by 


A 



dx dy {x, y) , 


where ca and Gy ^(x, y) are, respectively, discrete and continuous variational amplitudes, \AAJ’'T) 
are the NCSM eigenstates labeled by the set of quantum number A, ^y is an appropriate intercluster 
antisymmetrizer introduced to exactly preserve the Pauli exclusion principle, and 


\Ky) =[(|^ - «23 T’l) (|«2 a2/rT2)|fl3 

^ 5(x- 7723)5(i/-?7i,23) , 1 , 

^ 5 \ 1 / 

xm-i OT1,23 

are three-body cluster channels of total angular momentum J, parity n and isospin T where v 
represents a set of quantum numbers that describes the channel within the cluster basis. Here, 
\A - axi a\I’’^T\), \a2a2l2^T2) and |a3a3/j^r3) denote the microscopic (antisymmetric) wave func¬ 
tions of the three nuclear fragments calculated within the NCSM. The Jacobi coordinates describing 
the relative positions of the clusters are denoted by 7723 and 771,23- 

We calculate the unknowns of the NCSMC wave function [ca and Gy ^(x, y)] by solving the or- 
thogonalized coupled equations obtained by projecting the Schrodinger equation on the model space 
spanned by NCSM eigenstates and the NCSM/RGM basis. Those equations are solved by means of 
the microscopic R-matrix method in a Lagrange mesh |@] . Details on the procedure will be available 

in is. 

3 Application to ^He 

The lightest Borromean nucleus is ®He MM, formed by an "^He core and two halo neutrons. It is, 
therefore, an ideal first candidate to be studied within a three-body formalism. Hence, it was used as 
a test case when the NCSM/RGM formalism for three-cluster dynamics was introduced in 0 ] and 
here is studied again in order to perform a benchmark with such results. In this first calculation, we 







21st International Conference on Few-Body Problems in Physics 


0.05 

0 


-0.05 

q1 


-0.15 


5 10 15 20 

p(fm) 

Figure 1. Most relevant hyperradial contributions to the ®He g.s. wave function. Both the contribution from the 
NCSM wave function and the total NCSMC wave function are shown for a N„,„^ = 6 model space. The figure 
shows how the addition of the three-cluster basis within the NCSMC compensates the limitations of the NCSM 
to obtain an extended wave function characteristic of two-neutron halo nuclei. The hyperradial wave functions 
UfCyip) are the coefficients of the wave function when expanded in the hyperspherical basis, where K represents 
the hypermomentum. 
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describe the "^He core only by its g.s. wave function and couple the three-cluster basis with the ®He 
g.s. eigenstate obtained through the NCSM. 

We used the same potential as in lllll,i.e. , the similarity-renormalization-group (SRG) M 
evolved potential obtained from the chiral N^LO NN interaction 11 1511 with Agj^Q =1.5 fm '. With 
this soft potential the binding energy can be accurately computed by extrapolating the NCSM results 


to Nm 


2 , hence providing a good benchmark for the newly implemented NCSMC. 


From Table[T] we can see that the NCSMC g.s. energy quickly converges to the NCSM extrapo¬ 
lated value, unlike in the NCSM/RGM. This is due to the fact that the ®He NCSM eigenstate takes into 
account the short range correlations and "^He core polarization that are missing when considering the 
cluster basis alone. It is also important to note that, in contrast to the behavior offered by the NCSM, 
the NCSMC recovers the correct extended asymptotic behavior of the wave function. In Fig [1] such 
comparison is shown in a preliminary calculation at an Nmax = 6 model space. 

Finally, we can also compare the probability densities from the ®He g.s. obtained with the 
NCSM/RGM and the NCSMC. In Fig. |2] such comparison is shown and it is interesting to find that 
while the two main configurations (di-neutron and cigar) appear to have the same probability within 
the NCSM/RGM, the di-neutron probability is enhanced when using the NCSMC. This asymmetry in 
the strength of the probability peaks is known to be a characteristic of ®He and these results show that 
it is a consequence of the short range correlations. 


4 Conclusions 

The NCSMC uses an ansatz wave function that includes both an expansion in a continuum three- 
cluster basis and in a discrete basis of NCSM eigenstates. This provides a foundation that is capable 
of describing both short and long range characteristics of three-cluster systems. In the case of the 
®He g.s., we could see that this approach provides both the correct binding energy and extended 
asymptotic behavior unlike the NCSM that does provide the correct binding energy, but not the correct 
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Figure 2. Probability distribution of the ®He g.s. wave function in terms of the relative distance between the 
neutrons (r„„) and the distance between the center of mass of the neutrons and the ''He (ra,nn)- The di-neutron 
and cigar configurations appear to have the same probability within the NCSM/RGM (a), while the di-neutron 
probability is enhanced when using the NCSMC (b). 


asymptotics, or the NCSM/RGM that does the opposite. Calculations in larger model spaces for both 
g.s. and continuum states of ®He are underway. 
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